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20.  ABSTRACT 

The  hydraulic  system  frequency  response  (HSFR)  computer  program  was 
developed  to  simulate  the  dynamic  response  of  a hydraulic  system  to  the 
acoustic  noise  generated  by  the  pump.  Detailed  instructions  for  modeling  the 
system  pump,  lines,  and  components,  and  for  using  the  program  are  presented. 

For  a selected  system  pressure,  temperature,  flow,  and  pump  speed  range, 
the  program  calculates  the  pulsation  pressure  and  energy  levels  generated  by 
the  pump.  It  predicts  the  amplitude  and  location  of  the  resulting  acoustical 
standing  waves,  and  how  these  waves  are  transmitted  and  attenuated  throughout 
the  hydraulic  system.  The  program  may  be  used  for  acoustical  analysis  in  the 
pressure  side  or  both  the  pressure  and  return  sides  of  the  hydraulic  system..  s 
Estimated  line  lengths  and  sizes  from  preliminary  design  work  give  a 
good  estimate  of  hydraulic  system  natural  frequencies  and  pressure  amplitudes. 
The  program  outputs  plots  of  the  peak  flow,  pressure,  and  impedance 
amplitudes  of  any  selected  harmonic  of  the  pulsation  noise  versus  pump  speed 
for  selected  locations  in  the  system.  In  addition,  the  program  outputs 
plots  of  total  acoustic  energy  density  and  intensity  (power)  versus  pump 
speed . 
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SUMMARY 


The  hydraulic  system  frequency  response  (HSFR)  computer  program  was 
developed  to  simulate  the  dynamic  response  of  a hydraulic  system  to  the 
acoustic  noise  generated  by  the  pump.  Detailed  instructions  for  modeling  the 
system  pump,  lines,  and  components,  and  for  using  the  program  are  presented. 

For  a selected  system  pressure,  temperature,  flow  and  pump  speed  range, 
the  program  calculates  the  pulsation  pressure  and  energy  levels  generated  by 
the  pump.  It  predicts  the  amplitude  and  location  of  the  resulting  acoustical 
standing  waves,  and  how  these  waves  are  transmitted  and  attenuated  throughout 
the  hydraulic  system.  The  program  may  be  used  for  acoustical  analysis  in 
the  pressure  side  or  both  the  pressure  and  return  sides  of  the  hydraulic 
system. 

Estimated  line  lengths  and  sizes  from  preliminary  design  work  give  a 
good  estimate  of  hydraulic  system  natural  frequencies  and  pressure  amplitudes. 
The  program  outputs  plots  of  the  peak  flow,  pressure,  and  impedance  amplitudes 
of  any  selected  harmonic  of  the  pulsation  noise  versus  pump  speed  for 
selected  locations  in  the  system.  In  addition,  the  program  outputs  plots  of 
total  acoustic  energy  density  and  intensity  (power)  versus  pump  speed. 
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1.  INTRODUCTION 

The  Hydraulic  System  Frequency  Response  (HSFR)  program  predicts  how 

I 

oscillatory  flows  and  pressures  caused  by  the  acoustical  energy  content  of  a 
pump  output  are  transmitted  through  the  lines  and  components  of  a hydraulic 
system.  Resonance  can  occur  if  a frequency  of  this  oscillatory  output 
coincides  with  a natural  frequency  of  the  system,  or  of  any  part  of  the 
system.  Patterns  of  standing  waves  similar  to  those  observed  in  organ  pipes 
are  generated.  A resonant  condition  can  produce  large  oscillatory  pressure 
and  flow  amplitudes.  Resonant  acoustical  noise  can  cause  excessive  line 
motion  and  stresses  resulting  in  premature  failure  of  system  lines  or 
components . 

The  HSFR  program  predicts  the  pump  speeds  at  which  major  resonances 
occur,  and  defines  the  amplitude  and  location  of  the  oscillatory  pressure  and 
flow  standing  waves.  The  description  of  the  system  being  simulated  is  easily 
changed  to  investigate  various  practical  system  modifications  for  the 
attenuation  and/or  relocation  of  the  major  resonant  conditions.  This 
capability  allows  potential  problems  related  to  hydraulic  acoustic  energy  to 
be  eliminated  during  the  design  stage. 

The  program  is  written  in  Fortran  IV  language  for  use  on  the  CDC  6600  or 
IBM  360  computer  systems.  A punched  card  deck  is  used  to  input  the  program 
to  the  machine  while  output  is  by  line  printer. 

A user's  manual  for  the  HSFR  program  is  presented  herein,  describing 
in  detail  the  simulation  of  a system  and  the  interpretation  of  program 
output.  Volume  IV  of  this  report,  under  separate  cover,  contains  a detailed 
technical  description  of  the  HSFR  computer  program. 
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The  user  must  describe  the  system  to  be  simulated  by  means  of  punched 
data  cards.  The  description  includes  the  type  and  physical  characteristics 
of  each  of  the  elements  of  the  circuit.  An  element  may  be  a pump,  a section 
of  line,  a fitting,  a component,  or  a branch.  The  F-15  aircraft  number  one 
power  control  system  (PC-1)  is  used  herein  to  illustrate  circuit  modeling 
and  program  output. 

The  user  completes  the  problem  statement  by  specifying  the  range  of 
pump  speed  and  the  harmonic  of  interest,  the  locations  at  which  flow, 
pressure,  impedance,  and/or  energy  levels  are  to  be  plotted,  the  fluid  type, 
the  fluid  temperature,  and  the  steady  state  pump  output  pressure.  The 
program  calculates  the  applicable  fluid  properties  for  the  specified  system 
temperature  and  steady  state  pressure  via  a special  functional  subroutine, 
FLUID. 

The  program  calculates  the  oscillatory  pressures  and  flows  at  the  input 
to  system  elements.  Standing  wave  characteristics  produce  large  variations 
in  pressure  amplitude  along  the  length  of  a line.  Division  of  a length  of 
line  into  small  elements  may  be  required  to  allow  this  standing  wave  pattern 
to  be  examined  to  ensure  than  an  excessive  pressure  amplitude  is  not  being 
ignored 

Acoustic  analysis  may  be  performed  on  the  pressure  side  of  the  system, 
or  on  both  the  pressure  and  return  sides  of  the  system. 
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2.  PROGRAM  INPUT 


2 . 1 System  Description 

The  system  to  be  investigated  must  be  carefully  described  in  block 
diagram  form  before  the  data  input  cards  can  be  produced.  The  elements  which 
make  up  the  circuit  must  be  identified  by  type  and  located  sequentially. 

The  only  rule  which  must  be  applied  to  the  sequential  numbering  process  is  as 
follows : 


ALL  ELEMENT  WHICH  ARE  CONNECTED  IN  SERIES  MUST  BE  NUMBERED  SEQUENTIALLY, 
STARTING  WITH  THE  PUMP  AS  ELEMENT  NUMBER  1.  A BRANCH  LINE  MUST  ALSO  BE 
NUMBERED  SEQUENTIALLY,  WITH  THE  HIGHEST  NUMBER  APPLIED  TO  THE  TERMINATING 
ELEMENT. 


Figure  2-1A  illustrates  the  numbering  scheme  using  the  F-15  PC-1  pressure 
system  as  an  example  of  a complex  circuit  containing  several  branches  of 
various  lengths.  The  logic  of  the  program  would  deal  correctly  with  the 
same  system  if  the  numbering  scheme  in  Figure  2-1B  were  used. 

Figure  2-1C  illustrates  the  modeling  of  both  the  pressure  and  return 
sides  of  the  F-15  PC-1  system.  The  pressure  side  is  numbered  as  described 
above  up  through  the  terminating  pressure  element,  in  this  case  the  servovalve, 
element  #26.  The  next  sequential  element  (#27)  is  the  pump  inlet  which  is 
identified  as  a dummy  element,  NTYPE  7 and  KTYPE  1.  The  return  system  is 
then  modeled  sequentially  from  the  pump  inlet  up  through  the  terminating 
return  side  element.  In  this  example  the  return  side  is  modeled  back  to  a 
single  servovalve,  which  simulates  termination  of  the  return  side  at 
element  #36. 

Data  cards  must  be  stacked  in  order  by  the  user,  since  the  input  data 
does  not  include  the  sequential  element  identification  number.  The 
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FIGURE  2-1A 

F-15  PC-1  SYSTEM  DESCRIPTION  ELEMENT  NUMBERING  SCHEME 
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FIGURE  2-1 B 

F-15  PC-1  SYSTEM  DESCRIPTION  ALTERNATE  ELEMENT  NUMBERING  SCHEME 
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F-15  PC-1  PRESSURE  AND  RETURN  SYSTEM  ELEMENT  NUMBERING  SCHEMA 
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sequential  element  identification  number  is  printed  out  with  the  element 
input  data,  and  should  correspond  to  the  numbers  assigned  by  the  user  in 
describing  the  circuit.  The  program  is  presently  dimensioned  to  handle  a 
maximum  of  forty  (40)  circuit  element  records. 

2.1.1  Element  Description  - Figure  2-2  illustrates  the  element  descriptive 
data  for  the  illustrated  F-15  PC-1  pressure  system,  which  is  input  to  and 
printed  out  during  execution  of  a typical  program. 

Figure  2-2.1  shows  the  element  data  for  the  illustrated  F-15  PC-1  system 

combined  pressure  and  return  model.  All  the  types  of  elements  that  the 

program  is  capable  of  describing  are  discussed  in  Section  2.3,  including 

type  numbers  (NTYPE  and  KTYPE)  and  the  physical  characteristics  required  to 

define  each  type  of  element.  For  example,  from  Section  2.3.1  ELEMENT  NUMBER 

4 is  seen  to  be  a line  (NTYPE  1)  of  length  4.50  in.,  outside  diameter  0.625 

6 2 

in.,  wall  thickness  of  0.032  in.,  and  modulus  of  elasticity  16  x 16  Ibs/in. 
(titanium) . 

Each  element  data  list  has  space  for  a subtype  number  (KTYPE) . KTYPE 
allows  alternate  element  subtypes  to  be  described,  or  allows  subroutines  of 
varying  degrees  of  complexity  to  be  called  by  the  main  program.  Also,  KTYPE 
is  used  to  indicate  that  more  than  one  data  card  must  be  read  for  a giver, 
element.  In  the  example,  KTYPE=21  for  the  pump  (element  number  1)  indicates 
that  two  extra  data  cards  (total  of  3 cards)  must  be  read  for  the  pump 
(first  KTYPE  digit=2),  and  that  the  detailed  PUMP  subroutine  is  being  used 
(second  digit=l). 

2.1.2  Circuit  Termination  - Each  branch  circuit  must  have  a terminating 
element,  otherwise,  an  error  message  is  generated. 
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FIGURE  2-2 


Element  Input  Data 
For  F-15  PC-1  System 
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Elf »FNT 
*U**Pfc» 

**•• 

K 

PHYSICAL 

data******** 
data 

TYPE 

TYPE 

1 

0 

23 

.190 

.666 

1.120 

1.172 

,69« 

.67  0 

• 1*9 

.20000 

19,50000 

9.90000 

3.37500 

20,75000 

26.25000 

26.0 0*0  ft 

50,00000 

,06 ooo 

.76000 

2.20000 

.00042 

ftl),  ftpoon 

200.0-ftftftft 

2 

1 

-0 

6,000 

1.200 

.100 

30000000.000 

-o.ooo 

• ft  , ft  o 0 

- ft  , ft  ft  ft 

3 

1 

-0 

J.  500 

1.200 

,100 

30000000. 000 

-0,000 

• 0.000 
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<i 

1 

-0 

<Ooo 

,625 

.052 

16000000,000 

-0,000 

-0.006 

-ft, '*00 

5 

1 

-0 

3.  000 

,o2S 

,0.V 

16000000.000 

-0,000 

-o , rftft 

-ft.r  o 5 

6 

1 

• <> 

1 1 ,500 

.625 

,0  32 

1 6000000.000 

-c.OOO 

-0.006 

oft 

7 

1 

-0 

3,000 

.625 

,032 

16000000.000 
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5 

1 
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- '.ftftft 

1 1 

3 

-0 
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-0,000 

-0,000 

-0.000 

-0,000 

-ft  .ftftft 

1 2 

1 

•0 

6,00*’ 

.750 

,0  35 

lft'.COOOO.OOO 

- 0 , 0 c 0 

-O.oo ft 

- . 

1 * 

1 

• 0 

7 . 0 •»  0 

.756 

.0  35 

1 6 o 0 0 0 •)  0 . '*  0 0 

-0 . ftftft 

-' .000 

- . 

1 u 

6 

s 

•0,000 

-'*.000 

-0.000 

-0.000 

-0.000 

-r . ,»of 

IS 

9 

-0 

. 3 Uft 

-o.ooo 

-0,000 

-0.000 

— 0 , 0" 0 

- ft  . ft  ft  0 

lb 

1 

• l) 

2a, 000 

.500 

.029 

16000000.000 

-o.ooo 

-(>  # i.ftft 

. . ..4 

‘17 

1 

• ') 

369.000 

.375 

.01* 

16^000 00, 000 

- U , 0 o 0 

•6,000 

- . ft  A ft 

16 

1 

-0 

au , COO 

. 37S 

.032 

16000000.000 

-0.000 

-o  ,00ft 

- . ft  •!  0 

t« 

1 0 

-0 
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,7  70 

• o , 0 0 0 

-0.0%} 

-o.ooo 

- > ftf*  ft 

20 

u 

-0 

. 761 
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-0.0 00 

-ft.o 00 
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21 

1 

• 0 
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, 375 

.01* 
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-o.ooo 
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•ft ,f Oft 

22 

ft 

2 
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-c , ft  :• 

23 

1 
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2« 

19 

-0 
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2S 
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l« 
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-0.000 
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1 
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2* 
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-0 
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-0,000 
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• A ,00ft 
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FIGURE  2-2.1  Element  Input  Data  for  F-15 

PC-1  Pressure  and  Return  System 


A branch  circuit  may  be  terminated  by  any  element,  except  a pump.  A 
circuit  is  usually  terminated  with  a valve  or  a line  element.  Termination 


i 


with  any  element  other  than  a valve  results  in  a closed  circuit  analysis,  i.e. 
zero  output  flow.  A valve  termination  (NTYPE=14)  results  in  an  open  circuit 
analysis,  where  the  output  pressure  is  set  to  zero  to  simulate  an  open  ended 
termination . 

Circuit  termination  is  indicated  by  the  insertion  of  a "1"  in  front  of 
the  basic  NTYPE  number.  For  example,  element  number  23  is  a valve  terminating 
a branch  line,  consequently  its  NTYPE  identification  is  14,  rather  than  4. 

2 . 2 General  Control  Input  Data 

This  group  must  include  three  cards  to  input  the  title  of  the  run,  the 
number  of  elements  in  the  model,  the  fluid  type,  the  fluid  temperature,  the 
steady  state  system  pressure,  the  pump  speed  range,  and  the  harmonic  of 
interest . 

Card  1 - Title  Card 

Up  to  eighty  characters  may  be  used  for  a title,  which  will  appear  on 
each  page  of  output.  Use  a blank  card  if  there  is  no  title. 


K1S  PC-1  SYSTEM* PIJHP  Til  Rl.S  EREUUENCY  RESPONSE.  RLfi  M-'S  NllRMAI. 

| i i 1 | 5 b T » 9 i j il  | 3 I*  13  16  If  18  19  ?1  2l  Z2  75  24  | 2 |’  2B  :|:0  | 12  *3  3^9  M ST  l»  39  ft)  41  42  43  44  43  1 4 |d  43  *y  3 '.2  | VI  S3  !6  V ‘.<i  53  CJ  1.  6?  63  C4 
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} 1 3 3 ]'|]  J 3 3 3 3 ]|3j3|  3 3 3 3 3 1 3 3>3  1 3 3 1 3 3 3 3 3jl  1 3 1 3 ] 3 3 J 3'3  3 1 3 3 1 3 3 3 3;3  3 3 3 3 3 3 1 3 3;3  3 3 3 jjn  3 3 3;3  3 3 3 3 

1,1  I • | 1 ill 

4 4 4 I I'l  4 4 4 4|4  4 4 4 4 | 4 4 ||  I 4 4 4 1|4  4 4 4 4 4 4 4 |4'4  4 1 4 4 4 4 4 4 4'4  4 4 1 4 4 4 4 4 4 4 4 4 4 4 | 4 4 4 4,4  4 4 4 4 4 4 4 4 4 4 4 4 4 4 

1 y j 1 1 1 . t 1 y ig  :i  ti  11  1 ig  J'y  11  it  11  .■  .1  yy  )j  :i  jy|gy  ;i  11  ;y  ]*  ii  u ini  njis  ii  H n <t  41 1)  o 11  i.jil  hum  w .1  11  ;y  ,1  „[  i .r  « yt  it  ti  gy  ti  11 11 , it  1 II  II  1 >1  u 11  ,i  'I  >,  11  1 11 

3 i 1 5 I;i  i 5 i S i S 5 5 |;S  S 5 5 ! 3 S S i l|S  S S 5 5 5 1 S b |;|  i S S i | b S S |p  1 5 4 5 i 5 b S i|i  i 1 5 S S S S S SjS  5 S i 5 i S i 5 S’S  b 5 b S 

| 0 C E t<C  6 C 6 5 5 E 6 D C'C  i 6 8 6 G ( 6 | 6'S  6 C C 1 5 ( C C 6*  C 6 6 f 5 6 C |6'6  C 6 C 6 6 C C G C i C 1 1 C 5 C C G C'C  C G 6 C 6 6 6 5 G'6  G G 6 C 
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5 111  7 |;i  7 7 7 I 7 7 7 7 Tjl  ) 1 1 >1 1 7 J 7 ;?  7 7 1 1 7 7 7 ? 7]l  7 17  71)  ) 1 1 )jl  I 7 7 7 7 1 7 J i;i  7 I 1 I 7 7 7 7 7 ;7  71  7 I 7 7 7 7 7 |T  1 7 7 7 

2 I 18  8 l;a  ! 3 S 6 Hi  a i I '4 1 1 5 I II  8 I ijl  8 1 8 6 ! 8 1 8 l|l  3 1 8 n 8 8 ! M 8 1 8 6 3|S  I l;l  B 8 0 C I 8 8 I M 8 S I III  I IMS  8 81 

3 9 b 3 9)9  9 9 9 9 9 9 9 9 lib  9 9 3 3 9 S 9 9 9 1|  9 9 9 9 1 9 J 9 bib  9 1 9 1 3 1 9 9 j!l  9 9 9 1 9 9 9 9 b!l  3 !t|)l  9 3 l|l  9 9 9 9|J  9 9 9 9i9  9 9 9 9 


I J J • > ( 1 I 1 | I 1’  | TI  »*,  7,  1 1 ' ’»  ’1  ’|  i’  'I  '■  " 1 8 'I  ! 


« . »'  f (I  is  ft  1 1 14  H IJ  " " '1  I l ' I ’3  il 


Card  2 - Number  of  Elements,  Fluid,  Temperature,  and  Pressure 


Contains  the  total  number  of  elements  (NEL)  in  the  circuit,  including 
the  pump,  whether  the  pump  subroutine  or  the  empirical  pump  model  is  to  be 
used.  System  temperature  is  the  average  temperature  of  the  circuit  being 
analyzed.  Pressure  is  the  average  steady  state  pressure  in  the  circuit 
being  analyzed.  The  user  must  specify  the  type  of  fluid  used  in  the  modeled 
system.  A functional  subroutine,  FLUID,  described  in  Volume  IV,  computes  the 
required  fluid  properties  for  the  fluid  temperature  and  pressure  specified 
by  the  user.  The  program  computes  properties  for  MIL-H-5606B,  MIL-H-83282A, 
or  Skydrol  500B  fluids.  The  data  used  in  the  subroutine  FLUID  is  based  on 
the  best  source  data  available  on  fluid  properties.  The  user  may  also  input 
density,  adiabatic  bulk  modulus,  and  viscosity  data  directly  into  the 
program,  if  it  is  desired  to  use  data  or  a fluid  type  other  than  that 
provided  by  the  subroutine  FLUID.  User  specified  fluid  data  must  be  as 
corrected  to  the  specified  fluid  temperature  and  pressure.  Fluid  properties 
and  the  FLUID  subroutine  are  discussed  in  Volume  IV. 
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Card  2 - Continued 


COLUMNS  FORMAT 


DIMENSIONS 


11-20 


21-30 


31-40 


41-50 


51-60 


F10.0 


Total  number  of  circuit  elements  (NEL) 


User  option  = 0 


MIL-H-5606B=1 


MIL-H-83282=2 


SKYDROL  500B=3 


System  Temperature 


System  Pressure 


Fluid  Viscosity  or  blank 


Fluid  Density  or  blank 


Fluid  Adiabatic  Bulk  Modulus  or  blank 


IN^ /SEC 


LB-SEC2 /IN4 


(MIL-H-5606B  Fluiu  Specified) 
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■ 1 i i i lit 
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'i'i  i ' i i 
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Card  2 - Cont'd 


(User  Option-Direct  Fluid  Data  Input) 


30  0 100.  3000.  73.20H-03  34.00E-06  17.40E+04 
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' ‘J  <•  ' " ■'  "■*  " ? '1  '4  I'  '«»  If  H rj  20  n r.1  25  2«  25  2C  ?/  ;"?9  30  51  52  55  54  55  "•  if  50  59  40  41  42  45  44  45  40  4 1 40  49  30  51  W 55  54  35  5C  5*  *.  J V»  CJ 

oooq|qogoo<oo||0sooo  &;o  |||  o a c o 3 ofo  o|  o a o o o |fo  o |o  a o o o ofo  a|c  o o o o o;o  o c o o a o o o o.c  o o o o 
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Card  3 - Pump  Speed  and  Frequency  Data 
Contains  the  first  and  last  pump  speed  (rpm)  at  which  the  program 
calculations  are  to  be  performed,  the  increment  (rpm)  between  speed  calcula- 
tion  points,  the  harmonic  for  which  output  plots  are  desired,  and  the  number 
of  pumping  pistons. 

To  ensure  that  satisfactory  output  plots  are  obtained,  the  total  number 
of  increments  should  not  exceed  100,  the  start,  finish,  and  incremental  rpm 
values  should  all  be  even  numbers,  and  the  start  speed  should  not  be  zero. 
All  input  values  are  real  numbers.  If  no  harmonic  is  specified,  the  output 
plots  will  be  for  the  first  (fundamental)  harmonic. 


i 


« 
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COLUMNS  FORMAT 


F10.0 


DIMENSIONS 


21-30 


31-40 


F10.0 


F10.0 


Start  speed 


11-20  F10.0  Finish  speed 


Speed  increment 


First  Harmonic=l.,  0,  or  blank 


Second  harmonic=2. 


Third  harmonic=3. 


etc. =10.  maximum 


Number  of  pumping  pistons 


c-.'. ... 


. !.! 
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a 5 3 3 a so r 

1 1 ■. 
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2.3  Circuit  Element  Input  Data 


This  group  must  contain  one  card  for  each  element,  arranged  in  the  same 
sequence  as  the  elements  were  numbered  in  the  block  diagram  discussed  in  Section 
2.1,  The  pump  card(s)  is  placed  first,  followed  by  the  card(s)  for  element 
number  2,  and  then  by  the  remaining  cards  in  numerical  order,  up  through  the 
highest  numbered  terminating  element,  e.g.  valve  element  number  30  for  the  circuit 
illustrated  in  Section  2.1. 

The  data  required  for  all  of  the  types  of  elements  that  the  program  is 
capable  of  describing,  and  example  inputs  for  each,  are  discussed  in  the  following 
paragraphs . 
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2.3.1  Lines  hard-line  tubing 

. . . \ _ The  input  data  required  for  a 

„ o i i Ripid  Lines  (Tubing)  j.,i„D 

2. 3. 1.1  m n thickness,  and  Young  s modulus 

line  length,  outside  diameter, 

elercent  ate  line  same  lnp„t  data 

. ,r  , fitting  element  may  be  modeled  usi  8 
of  elasticity.  A titting 

parameters  as  for  a rigid  line. 


' 


2-15 


2. 3. 1.1  Rigid  Lines  (Tubing)  - Cont'd 


COLUMNS 

FORMAT 

1-5 

15 

6-10 

15 

11-20 

F10.0 

21-30 

F10.0 

31-40 

F10.0 

41-50 

F10.0 

DATA 


NTYPE  = 


KTYPE  = 0 or  blank 


Line  Length 


Outside  Diameter 


Wall  Thickness 


Modulus  of  Elasticity 


DIMENSIONS 


I 23456789  iQi 


. Odd  160001100. 

J li  U O )<  S3  SC  Sr  IS  39  43  41 


■4  .3  16  .7  10  13  20  2#2  23  .4  23  2C  2 7 28  r.  fl  li  U 43  S4  S3  SC  37  IS  39  43  4i  44  4J  «g  47  43  44  33  31  32  33  14  35  34  57  38  39  CO  61  02  63  t4 


. 2 S 4 3 6 7 * *,  10  II  12  13  14  13  tO  17  10  19  20  2i  2?  23  24  23  26  27  28  29  SO  Si  3?  33  34  S3  SC  37  38  39  40  4i  42  43  44  43  40  47  40  49  50  31  32  53  54  ",  V.  • • 3»  **»  6f 

0 0 8 0 OiO  0 0 0 0 0 0 0 0 0-0  0 0 0 0 0 0 0 0 O'O  0 0 0 0 0 0 0 0 O'O  I 0 0 0 0 0 0 linn  o q a 0 0 o;o  0 a a 0 3 0 0 0 o;o  o o o c a o o o o;o  q o o d 

1 1 1 1 iH  ! 1 t 1;  1 11  (i  11  : ' 1 1 1)  it  it  it  ji  :i  n tt  n|tt  ti  11  it  It  it  :i  u :i  njjt  tt  it  11  n 11  *t  «J  it  ti.it  it  11  it  it  11  it  it  it  n|it  it  i.  it  w ti  it  tt  ti  t'.|it  it  ti  t,  tt  n it  it  11  lift  tt  n . it 

I I I I |'  1 I I I I I I II  1 <1  II  II  I ! ! Il'll  I I II  I II  Ijl  I I 11  I | 11  Ijl  I 1 I I I I I I 111  1 I I 1 I II  PI  I I I I I I II  Ijll  I 1 I 

2227  2<2  22222222  2l2  22222222  2!2  2|222222  2>2  22|22222  2,2  22222222  2<2  22222222  2'2  22222222  2,2  2222 

1 , 1 “ 1 1 1 1 1 

1 , 1 1 1 1 1 1 

1 3 3 2 Id  3 3 3 3 1 3 1 3 3fl  3 3 3 3 i 3 3 3 3 1 3 3 3 3 1 3 3 3 3|  ]|  3 3 3 3 3 3 3|3  3 3 3 1 3 3 3 3 3;3  3 3 3 3 1 3 3 3 3;3  3 3 3 3 3 3 3 3 3;3  3 3 3 3 

4 4 4 4 4<4  44444444  |'4  4444  3444  4<4  4444  4 444  4*4  44444444  4*4  44444444  4'4  4 4 4 4 4 4 4)4  44444444  4 M 4 4 4 4 

I | J I i|l  I 1 1 1:  II  * ) il  II  ,|i»  '!  I U !rl  72  75  M J'j J’l  21  21  71  K 1'  12  1)  14  )»'(,|  )l  II  11  48  (l  :}  4)  44  «/(4S  41  4J  41  St  .1  1?  •.]  U JS,K  W 31  l.  ii  12  II  14  tSi»;  «l  (I  IS  0 H 12  I)  14  \ v 1 I r|  || 

5 5 5 5 5'5  5 5 5 5 5 5 5 5 Ir3 1 5 5 5*4  b 5 5 5'5  5 9 1 b 5 5 5 5 5<5  5 5 5 5 !.  1 5 5 515' 5 5 5 5 5 5 5 5 5'5  5 5 5 5 b b 5 5 i'5  b 5 5 5 i 5 5 5 5‘5  5 5 5 J 

I 'll  I ' 1 1 

I « I I « I * * ' 

t b 6 6 E 6 6 6 6 $ C G 6 G C'C  C 6 6 616  C C C 6'G  1 1 C C 6 C C C 61C  G G 6 6 ' J 1 6 U 6 6 G C 6 C G C t'C  6 C C C C C C 6 S C 6 C 6 6 0 C 6 t 6^6  8 6 6 C 
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1 1 1 n;i  1 n 1 1 1 n I'm  i n n ii'M  m 1 n 1 I'm  n 1 ; 1 11 1'l  1 n 1 1 m i;i  n 1 1 m ii'i  m 1 1 1 1 1 ti ) 1 1 1 

1 1 1 1 1 ' 

8 8 8 3 8'8  8 8 8 8 ; 9 8 8 !|  3 0 3 8 8 8 8 3 8 | 0 8 8 8 ! B 8 I l'|  6 9 8 8 1 3 8 8 8'!  8 8 8 1 8 8 8 8 B’8  8 8 8 I 8 8 8 I SI  8 8 8 8'8  8 8 8 I'S  8 B B 8 

• j 1 ' l I • III 

it  1 1 i!i  n 1 Mi  i h s!i  1 s h 9 m n!n  s 9 1 n m!i  1 h s !.  1 m!i  n 9 > n n i!i  9 1 9 9 m 1 t!i  1 1 n!i  n n!n  m 


2. 3. 1.2  Flexible  Hose  - The  input  data  for  a hose  are  the  active  length  excluding 
the  hose  end  fittings,  the  inside  diameter,  and  the  volumetric  expansion  due  to 
pressure  expressed  as  an  equivalent  bulk  modulus.  Equivalent  bulk  modulus  is 
obtained  from  experimental  data  centered  about  the  working  pressure  according 
to  the  following  expression. 

Bulk  Modulus  = (Pressure  Change)  x (Total  Hose  Volume) 

(Volume  Change) 


columns!  format 


15 


6-10 


11-20 


21-30 


31-40 


41-50 


15 


F10.0 


F10.0 


F10.0 


F10.0 


DATA 


NTYPE  = 1 


KTYPE  = : 


Line  Length 


Inside  Diameter 


Blank 


Bulk  Modulus 


DIMENSIONS 


1 1 43.0 

J « a i T * 5 tO  ii  12 1 -4 


.7*50  1 00000. 

15  15  i T 18  15  25  2S  S4  25  2&  IT  21  25  10  Si  U 55  14  15  5.  IT  58  55  40  41  44  45  45  4T  40  41  55  5l  5?  SS  54  55  55  5T  54  55  US  t ■ 52  5 S 


9 « 8 fc  r * 'i  10  Ii  12  ij  14  15  iv  i f 18  19  20  2 1 Zi  25  24  25  26  »7  28  29  50  5i  52  55  34  S3  3t  37  38  59  40  41  4i  43  44  4*.  4b  47  40  «9  30  ii  i.  35  34  33  •* 
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Ji 

I 

' 

I 


2.3.2  Resonators 

2. 3. 2.1  Lumped  Volume  Resonators 

The  input  data  for  a lumped  volume  resonator  includes  the  length  of  the 
neck,  inside  diameter  of  the  neck,  neck  wall  thickness.  Young's  modulus  of 
elasticity  of  the  neck  material,  and  cavity  volume.  Entrance  angle  effects 
are  not  considered. 


2-19 


2. 3. 2.1  Lumped  Volume  Resonators  - Cont  d 


6-10 


11-20 


21-30 


31-40 


41-50 


51-60 


FORMAT 

DATA 

15 

NTYPE  = 2 

15 

KTYPE  = 0 or  blank 

F10.0 

Length  of  Neck 

F10.0 

Internal  Radius  of  Neck 

F10.0 

Wall  Thickness  of  Neck 

F10.0 

Modulus  of  Elasticity 

F10.0 

Cavity  Volume 

DIMENSIONS 


£ 3.0  1.000  .100  30000000.  13. 

i ? : « i «.  r h t it>  1 1 i 1 14  i5  i«  ir  it  io  zo : |;  ?5 2«  25  :c  sr  :•  ;i  | n u ;i  i«  n jc  57  :o  S3  40  41  -i|  ; 44 « 40  «.-li  49  50  n 52  11  44  44  h ir  jisim  6>  «2  u *« 


lb  17  10  19  20  2l  22  25  24  24  26  it  JC  29  50  51  52  55  54  55  56  5 2 58  59  40  4i  42  45  44  45  44.  4F  40  4'i  tO  51  V»  55  54  54  . 


0 0 0 0 0;3  3 0 3 010  0 0 0 GO  |0  Q C 3 0 0 0 0'3  in  0 0 0 0 0 0|3  0 II  GO  0 mini  o p o o o g;o  o o o o o o o o o;o  o o o o o a o o c;o  o c o 

1 i i i 5*i  i i i cl  i ■,  ; ii  nj  i ii  i ii  :i  n ??  jj  ?i  »]it  v ?i  ?»  n n v » u j/,;*  u n n u h o «i  <«  n * » :t  51  y.  '4)  51  .sJm  -i  j \i  >o  11  kj  ii  i<  11,(1 11  0 u 1 ,1  11  n u «vi  11  1 n u 

1 I II  l|l  II  II  I I I I Ijl  1 I 1 I I 11  I 1 1 II  I I I I I I Ijl  | I I I I 11  I ijl  11  I I I II  | Ijl  I I I I I I II  Ijl  I II  II  I I li; I I II  I 

2 2 2 2 f?  7 2 7 7 2 7 2 2 2^2  2 2 7 2 2 2 2 2 2l2  2 2 2 2 ’ 2 2 2 2^2  2 2 2 2 2 2 2 2 2'2  2 2 2 2 2 2 2 2 2^  2 2 2 2 7 2 2 7 2'2  2 2 2 2 2 2 2 7 2'2  2 2 2 2 

3 3 3 3 3 J 1 3 3 3 3 3 3 3 |'|  3 3 3 3 J 3 3 3 3l  3 3 3 3 3 3 3 3 3 | 3 3 3 3 3 | 3 3 3'3  3 3 3 | 3 3 3 3 :-'|  3 3 3 3 3 3 3 3 30  3 3 3 3 3 3 3 3 30  3 3 3 3 


4 4 4 4 4 I 4 4 4 4 4 4 4 4 4 '(  4 4 4 4 4 4 4 4 4'4  4 4 4 4 ( 4 4 4 4 ' 

1 * I 1 . j<  1 1 t k'  1 w ii  ii  n|  i ii  11  11  ii  11  :>  ;i  1 :ij;i  :i  ;i  n i:  n 1?  n )i  »] 

5 5 5 1 io  5 5 5 5 1 5 5 5 5 '5  5 5 5 5 5 5 5 5 5 S 5 5 5 5 5 5 5 5 5' 

i'll 
I I I I 

e 6 6 6 6 3 6 6 6 G G C G 6 G 5 6 6 G G 6 S 6 G 6<5  6 C 6 G 5 G G C S' 
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I 
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i 
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2 J I 5 I ' I I 'I  ll  ? II  il  '5 
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2. 3. 2. 2 Pulsco  Acoustic  Filter 

The  input  data  for  an  acoustic  filter  is  on  four  cards.  The  first  card 
contains  the  three  volumes.  Each  of  the  last  three  cards  contains  the  length, 
outside  diameter,  wall  thickness  and  line  modulus  of  elasticity  for  each  of 
the  three  internal  lines.  Line  data  is  input  in  the  order  of  line  1,  line  2, 
and  line  3. 


NOTE 


This  device  is  manufactured  and  marketed  by  the 
Pulsco  Division 

American  Air  Filter  Company,  Inc. 
Louisville,  Kentucky 

The  design  and/or  inventions  disclosed  are  the  property 
of  American  Air  Filter  Company,  Inc.,  Pulsco  Division. 


2. 3. 2. 2 Pulsco  Acoustic  Filter  - Cont'd 


First  Card 


6-10 


11-20 


21-30 


31-40 


FORMAT 

DATA 

15 

NTYPE  = 2 

15 

KTYPE  =32 

F10.0 

Volume  A (VA) 

F10.0 

Volume  B (VB) 

F10.0 

Volume  C (VC) 

DIMENSIONS 


• 0 0 0 0 .8  0 !l  0 0 

0 0 0 C C '0  0 I !l  1 0 3 0 0 0 

■0  0 1 0 | 
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2. 3. 2. 2 Pulsco  Acoustic  Filter  - Cont'd 


Cards  2,  3,  and  4 (Typical) 


COLUMNS 

FORMAT 

DATA 

DIMENSIONS 

1-10 

F10. 0 

Line  Length 

IN 
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IN 
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2.3.3  Volumes 

2. 3. 3.1  Volume  Without  Losses 

Input  data  for  a volume  element  is  the  volume  itself.  This  assumes  no 
internal  losses  in  the  volume  element.  For  example,  a filter  cavity  may  be 
modeled  as  a simple  volume. 
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2 .3. 3.1  Volume  Without  Losses  - Cont'd 
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11-20 
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NTYPE  = 3 
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2.3.4  Valves 

2. 3.4.1  Simple  Valves  - The  input  data  for  a valve  element  includes  the  valve 
gain  linearized  at  the  steady  state  circuit  flow  through  the  valve.  Determination 
of  valve  gain  is  discussed  in  section  2. 3. 4. 2.  Valve  gain  is  expressed  as  pressure 
drop  (psi)  per  unit  flow  rate  (cubic  inches  per  second) . Flow  out  of  the  circuit 
being  analyzed  is  input  for  terminating  valve  elements.  This  "overboard  flow" 

is  the  steady  state  flow  through  the  terminating  valve  at  the  input  steady  state 
system  pressure. 

2. 3. 4. 2 Determination  of  Valve  Gain  - Typical  valves  which  may  be  modeled  are 
system  shut-off  valves  (low  gain),  electrohydraulic  servovalves,  mechanical  servo- 
valves, and  combinations  of  electro-mechanical  servovalves,  such  as  may  be  found 
in  an  integrated  actuator  package. 

For  a valve  pressure/flow  relationship  of  the  form 
P = KQn 


where:  P = pressure  drop  (psi) 

Q = flow  rate  (cis) 

K = constant 
n = flow  exponent 


the  linearized  valve  gain  (G)  may  be  determined  from 


G = 


dP 

dQ 


nKQ: 


,n-l 


but  K = ^ 

, c „ nP 
therefore  G = — . 


P = KQ 
dP 


dQ 


= G 


Steady  Stat^ 
Point 


If  the  valve  flow  can  be  characterized  as  an  orifice  (n  = 2) , then  the 
2P 

gain  is  G = / . The  orifice  relationship  is  typical  of  electrohydraulic 

servo  valve  steady  state  control  flow.  If  the  valve  flow  can  be  characterized  as 

p 

laminar  for  the  steady  state  condition,  then  n = 1 and  the  gain  is  G = / . The 

laminar  relationship  is  typical  for  null  leakage  flow  across  lapped  spool  valves, 
e.g.  mechanical  servovalves,  and  the  second  stage  of  an  electrohydraulic  valve. 
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Parallel  valve  elements,  for  instance  those  within  an  electro-mechanical 
integrated  servoactuator , may  be  combined  for  modeling  as  a single  valve  element 
by  computing  an  equivalent  gain  (Ge)  for  all  the  parallel  flow  paths. 

-i-  *-  + *-  + . . . 

Ge  G1  G2  G3 

Empirical  pressure  drop/flow  data,  if  available,  should  be  used  to  calculate 
the  gain  at  the  steady  state  flow  condition.  The  flow  relationship  may  be 
assumed  unless  the  flow  exponent  is  available  from  empirical  data. 
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2.3.4  Valves  - Cont'd 


COLUMNS  FORMAT 


DIMENSIONS 


11-20 
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2.3.5  Accumulators 

2. 3. 5.1  Piston  Type  Accumulators 

The  input  data  for  a piston  type  accumulator  are  the  piston  mass, 
cylinder  internal  radius,  cylinder  wall  thickness,  cylinder  modulus  of 
elasticity,  maximum  gas  volume  (piston  bottomed  to  oil  side),  minimum  gas 
volume  (piston  bottomed  to  gas  side) , gas  precharge  pressure  (piston  bottomed 
to  oil  side),  working  pressure,  and  viscous  friction  damping.  Two  data  cards 
are  required. 

Viscous  friction  damping  data  is  not  readily  available.  Its  value  will 
be  determined  experimentally;  in  the  meantime  an  estimated  value  will  be  used. 

• 


i 


i 
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2. 3. 5.1  Piston  Type  Accumulators  - Continued 
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1 8 C 8 C 3 S 8 8 8 

| 8 8 d 8 

9 8 8 8 7 | i,  8 8 8 o ! 

5 ! J ! 8 8 J 9 9 9,9  8 ! 3 9 J 9 9 3 9 

3 9 3 9 0 3 9 9 9 

99995  988999898 19  3889  ;9  9999 

9993999999 

• a s 

39539  !9  1999. 9 3999  5 

, 1 1 1 1 t 1 » 1 « '» 

II  1 1 IJ  <4  Ii  H H 1 '1  31 

fi  n i ’i  ’i  i 

• )’  J)  A 2 It  !'  It  • * I * 4)  ll  43  4«  i ,1  »•  M - M V 

r>  «>  -, 

I 
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Second  Card 


F10.0 


F10.0 


Working  Pressure 


Viscous  Friction  Damping 


LBS ' SEC/IN 


3 0 0 0 . 


9 * 1 'i i i it  It  as  21  !0  SI  5?  53  5«  J5  36  17  it  39  «0  4i  42  43  44  45  46  4?  40  49  60  5:  52  53  54  55  5C  5?  50  59  Cn  6'  6?  61  04 


• 2 1 4 5 c.  I h 9 i„i  12  1 1 4 i !,  ■ 0 i 1 iu  >J  20  ? 77  ! :»  ?5  ?6  /l  2 ) 23  JO  Jl  J?  S3  34  33  36  37  38  30  4 0 4,  47  4*  <4  45  15  4/  40  41  50  5i  5?  *3  54  V.  *6  *7  *3  59  60 

0 ii  1 1 1 n c o c 5 o o o | j o o o o a o o o Q'3  o ii  o o o o o o co  o o o o o o o c o;o  o c o c o ii  o o o;3  o a o o a o g o o;o  o o c oj:  o o o o;o  o o o o 

'll*  ' i ii*  ( i ii  ii  n h i j;  ji  hjis  it  ;i  n jo  .i  j?  i j ,«  uju  jj  :i  n :»  n u u u <i,«t if  u is  S9  n sj  si  .1  ssjn  .1  ,1  os  u 11  v u n 11,11  v 1:  ii  ic!  1 12  n 11  »Si  1 m i s ij 

1 I II  l|lll  III  I 1 1 H | H I I i ! 1 M III  II  I I I I'l  III  1 III  Il!l  111  III  fll'I  I I I I III  M'l  II  I ill  I lll'l  I I I I 

. 1 1 j 1 'I 

2 2 2 2 2*2  72227722  2*2  22227212  2*2  22227227  7*2  22222222  2'2  22222222  2‘2  22222222  !*2  2 2 2 212  2 2 2 2‘I  2 2 2 1 


2 2 7 2 2'2  2 2 7 21?  7 2 7 7 


]| n )| 3 } 3 ]| n n 3 

I | 

444414444414444 

I 2 : i i{l  i I 9 J.  ; « 

5 5 5 5 5*5  5 rj  5 S tb  5 5 b b 

1 

6 t G E li  6 G C 6 'i  6 6 (,  i 

: I 

j 1 1 1 i;m  1 1 1 1 1 j i 

• I!  3 S I I a 8 J i ■ i 8 1 


unuiinunii 


i ' i ii  20  > :/  ;;  u j;(;s  <>  ;i  n jc  ,•  j?  ::  ,t  j> 

1 1 1 1 1 1 1 1 1 m 1 1 1 ii 1 1 1 1 

i 

272227722  2>2  222222277 


j|  3 3 3 3 3 3 3 3 3|3  3 3 3 3 3 3 3 3 3 

444444444  4'4  444  44  1444 

|'i  H » ii  ii  it » n :»jii  ii  ,i  ;i  *jii  ii  u » ,, 

b i b b b i i b b blb  b i i i i i b b b 

i 

.ubGtjCCE  65  CC6L6G66E 

M CDONNCLL.  AUTOMATION 

;i  1 1 ; 1 1 1 ill')  1 1 m m i 

I B 4 Hi  i U II  Ii  I)  0 8 a 8 8 8 8 3 


1 3 3 3 3 3 3 3 3 3'3  3 3 3 3 3 3 3 3 3 

i 

444444444  4*4  444444444 

.6  H it  JO  (3  II  4?  Ii  II  4. Jit  II  l|  II  31  |J  34  0) 

5b!)bS!)S‘  55:5b55S'j555S 
\ 

: 6 6 C G b G ' l G*C66CGS666C 
COMPANY  [ 

1 1 11  71?  ? : MJI  I 1 1 7 1 1 7 1 1 


im  oh  0 'J  9 'J|*J  9 9 9 9 |'J  9 9 9 9,9  9 9 9 919  1 . '•  3»9  9 3 9 9|3  9 9 9 9 


31  ’ l A IS  Ii  II  V 11  li  tb.tl  I IS  tl  ic!  I !J  I)  II  ISl 

I l l I I I l l l I ' I I l I l - [ M I I 
: | : 
272222222  2'2  2 2 2 2 2 2 2 2 21 


3 3 3 3 3 3 3 3 3 3'3  3 3 3 3 3 3 3 3 3'3  3 3 3 3 


444444444  4 14  4 4 4 4 4 4 4 4 4 

'i  if  fl  1}  l.  t tf  li  Ii  t. ,it  I;  II  [1  ,1  ll  II  ,1  u i. 

b 5 5 5 5 1 5 b b bjb  b 5 1 5 b b b b b 

' 

0 C G 6 5 6 C C G GjG  6 0 C 6 |b  6 6 G E 

1 I 1 I I I I I I l;l  I I I I I I i 1 I 


I 1 8 8 [ 8 0 . i S 8 1 8 8 8 9 8 8 1 0 'S  8 8 8 Bid  S 0 I H 0 8 8 S 8 3 I 8 I 


3 9 9 9 3|J  9 9 9 9,9  9 9 9 9'3  9 3 3 9,9  9 9 9 9 


i i i . i ,ii 


i ii  t -i  ii,  ' i 
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2.3.6  Branch  Element 


The  data  input  required  to  describe  a branch  element  is  the  total 
number  of  elements  in  the  branch.  For  instance,  in  Figure  2-2  element  number 
18  is  a branch  with  5 elements  in  the  branch,  element  number  26  is  a branch 
with  2 elements  in  the  branch.  The  number  of  elements  in  a branch  is  used 
as  an  integer  number  by  the  program. 
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2.3.6  Branch  Element  - Continued 


FORMAT 

DATA 

15 

NTYPE=6 

15 

KTYPE=Number  of  Branch  Elements 

DIMENSIONS 


i i a i1)  ;o  it  is  it,  n :•  :*  so  31  si  is  si  is  sc  :r  m s»  40  41  4!  «s  77  15  «t  47  41 79  so  n s?  ss  si  ts  s»  sr  si  ss  to  n t;  u .« 


‘Ill) 

Mill 
7 177  7 
11  3 3 3 

4 4 4 4 4 

1 1 1 1 i 

5 5 5 5 5 

6 6 C C | 

mu 

h m 


Ji  ?»>  ?7  29  JO  3'  3.’  JJ  34  33  Jb  J.>  ja  39  40  41  47  43  44  43  4fc  47  49  49  J.0  3i  W ?J  34  33  ! 6 V •_*  »•}  t.  j 

0 0 q 0 o a 0 q 0 olo  0 0 0 a a 0 o q c|3  0 a 0 Qio  0 0 0 o!o  a o a o;o  0 Q 0 0 3 0 0 0 c;o  0 0 0 o 0 0 0 0 co  0 0 o 0 0 0 0 0 o;o  0 0 0 0 0 0 0 0 o;c  a o o 0 

? 7 ( 1H  I I I ic!  7 I n ll!'i  II  7 IS  ?t|7l  II  7,  71  :i!ll  II  II  71  il|  I 77  77  II  ill  II  Jl  71  U II  77  77  II  IS.II  II  II  II  SI  SI  17  SI  St  Ssjll  SI  SI  SI  H II  II  17  M IS,«  II  Is  II  l|  II  l|  II  ;|  isi’l  I I 77 

I I pi  l 1 t I l 1 l 1 l!l  t 1 1 1 1 1 1 1 ft  I I I I I I I i l'l  l l l i l 1 1 i i'i  l i i i 


i 77  II  77  7 i 7 )'  il  77  'C  II  1'  7]  77  71.11  1 1 II  77  SI  SI  17  SI  SI  SS  )i  SI  SI  SI  11  II  17  I]  17  is,ll  II  Ij  cl  ill  ' 17  17  I ' 


I I 1 I 1 I I I 1 M I I I 1 I M I l'l  I 1 I I I I II  l'l  I 1 I I II  I I I I I I I I I I I I M I I 11  I I I 1 l'l  I I I 1 1 I I I 1 

i i i i ■ i 

222227222  2 12  32222222  2*2  22222722  2‘?  22222222  2*2  22222222  ?'2  22222222  2'2  222222222 

, i i i i i 

m n i n : in  1 1 1 1 1 n i'i  1 1 1 1 1 n i i'n  n i n i n'i  1 1 n n n i'i  ni  1 1 3 n i'i  3 n i n n i 

i i i i i i 

i i i i i i 

4 4 4 4 4 I 4 4 4 4 i 4 4 4 4 I 4 4 4 4 4 4 4 4 4 4 4 4 4 4'4  4 4 4 4 4 4 4 4 4 M 4 4 4 4 4 4 4 4 4'4  4 4 4 4 4 4 4 4 4'4  4 4 4 4 I 4 4 4 4 

71  ' '7  II  '1  IS  |‘l  I '7  II  It  .'7  71  *i  Isjll  I *1  'I  17  II  17  II  77  nil.  7' II  77  I*  71  1?  «3  11  7l[  . 71  71  17  St  SI  II  S7  SI  isjsl  SI  si  si  II  .1  71  II  17  islci  il  H 71  11  M 17  II  ;l  ' 

5 5 5 5 ||5  5 5 5 5,5  5 5 5 5 5 5 5 5 5|5  5 5 5 5 5 5 5 5 5|5  5 5 5 5 5 5 5 5 5|5  5 5 5 5 5 5 5 5 S|5  5 5 5 5 5 5 5 5 5|5  5 5 5 5 3 5 5 5 5 

‘ ' ! ! 1 1 I 

6 5 6 G 6 0 S 5 0 6 6 6 5 G C 3 0 C C C>C  6 5 0 6 16  6 G 6 6 1£  5 S 6 C 6 6 r,  6 6 >a  C C 6 C 3 G 6 E £'6  5 0 C E 6 6 C 6 U 5 C 0 C 6 5 6 G C 

MCCJONNELL  AUTOMATION  COMPANY  J | 

1 1 1 1 1 m i M 1 1 ; ill  1 1 iM  1 1 1 ill  1 1 1 i'i  1 1 1 1 1 1 mw  1 1 1 1 1 1 1 1 ri  1 1 1 1 1 1 1 ? i'i  1 1 1 1 ? 1 1 1 1 

iii  i < i i 

3 8 3 8 13  ! 6 8 8 3 S 58  8 3 8 0 8 IS  I 3 8 6 J 8 8 B ! 'j  S3  3 £ 8 l 3 5 8 'i  i S I S 3 3 ! J I B 8 ! 6 6 3 II  8 I 8 's  8 8 8 1 J 8 8 9 8 


7 7 7 1 7 

mu 

4 4 4 4 4 

<i  n ii  ii  ii 

5 5 5 5 5 

8 6 6 C B 

1 1 1 1 1 ' : 


n 3 9 Sl9  1 9 9 9 9 3 1 3 5|9  9 9 3 9,3  9 9 9 013  9 i 9 9,9  0 5 ! 5|J  ! 9 9 Jit  5 J J 9|9  9 9 9 9,9  3 9 9 S|3  9 9 9 9,3  9 9 9 9|9  9 3 3 S. 


1 I S I ' I I I II  ll  '1  17  Is  1 ' I 71  71  7’  II  II  11  ’i  'I  'I  ’7  II  1 71  77  IS  1.  I 7:  'I  7*  ' •'  ii  ■■  IS  " 1 


>l  11  «•  n »,  M (»  (I  «»  H • 7} 
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2.3.7  Dummy  Element 

Input  data  required  to  represent  a dummy  element  is  NTYPE=7.  Capability 
for  the  use  of  a dummy  element  is  sometimes  convenient  to  minimize  element 
numbering  changes  for  describing  similar  circuits.  Subtype  1 is  used  to 
designate  a pump  inlet  when  modeling  the  return  side  of  a system. 


MAC  2S369  (2’  *UL  39> 


2.3.7  Dummy  Element  - Continued 


COLUMNS  FORMAT 


DATA 


NTYPE=7 


Pump  i n 1 e t KTYPE=1 


All  others  KTYPE=0  or  blank 


DIMENSIONS 


Pump  Inlet 


p i ',)Tf;am  c>  aif.mfn: 

I i<  *■ 

0!»  0 0 P I 1 0 : POOL'  C ,1  0 0 0 0 P (! !!  Ii  u 3 1 0 fl  C C 0 fl  1 D B 0 0 0 : i)  0 (1 i 0 n 0 0 0 " !)  3 » ? ' 0 f " 3 
•Ij  ■ | : ” * »•»•*••  ’ i •'  ••••■•  ' ' 

in  i i iii  i i i|  1 1 i i i i i i I i i i t i i i i i i i i ; i i i i i i i 1 i i ' i i i i m i i ' i i i i * 1 


III  1 

\ 

1 i ! 

712  2 
| 

? 7 ; 

ill  1 

3 3 i 

1 

il;  i 

4 

1 

1 

HO  '• 

1 

r,  ii' 

1 

SIR  0 
1 

ij  C G 

/I/  1 III 
1 

1 

Cl ^ 5 

3 8 4 

ill  0 
V 1 

38* 

* 'l 

'7  7 2'.*: 


i i 1 i 1 I i i . ; i 

7 7 7 ' ; 7 V * 7 


1 I in  II  1 I 

i 

7 i 

i 


.in. 

.'  0 5 5 


• • i 1 3 1 .•  3 3 n : 3 3 3 n 3 2 j : -i  .y  * i *• 

*!  1 ; I 1 4 1 4 4 4 ' 4 •'  4 4 4 • 4 ; 14  4't  4 4 4 


iii  r /in  i n m 1 1 1 n i n 1 1 1 m m 1 1 1 1 n 1 1 m 


r v.  6 «.  e f r r-  s o f,  & 6 : 

i n i ' n n n * 


i 8 3 9 0 5 8 3 8 9 3 V:  i ' * 3 : 
i r a r,  (•  s 5 c g s tin  mi- 


i ; 1 1 ii  l n i ir  1 1 1 1 


318  8 3 8 > 8 3 8 8 8 8 8 8 3 8 8 8 8 M 8 .U  4 8 8 8 8 I 8 0 8 8 8 8 8 8 8 8 e 8 3 8 ? P 8 8 r 8 5 S 0 3 8 8 8 8 - a s 5 8 5 •.  4 8 8 

ijg  035  3010391339  !)99*1''99999‘)(tli'll!19973'!!19D  '4  t t o n ^ n • ‘i  3 9 3 VI  i ' ' ’ ' 7 3 C 1 a 'l  2 

• I?  j i • . I > « *•  -mi  »j  i»  ii  . imj  ii  * ’ * ’ * * 1 " • i «•  i* « •:  • ’ • * 


. 9 liK  3 " V 


2.3.8  Quincke  Tube  - Input  data  for  a wide  band,  helical  Quincke  tube 
resonant  mode  damper  are  outer  (delay)  tube  length,  inner  tube  inside  and 
outside  diameter,  outer  tube  inside  and  outside  diameter,  wound  element 
cross-sectional  area,  helix  pitch,  number  of  holes,  hole  length,  hole 

spacing  distances,  and  hole  diameters.  The  maximum  number  of  holes  is 
16  with  the  existing  program. 

This  device  is  shown  schematically  below.  The  outer  spiraled  passage 
for  the  illustrated  model  is  formed  by  winding  a solid  element  around  the 
straight  inner  tube,  which  is  then  enclosed  in  another  straight  outer  tube. 


2.3.8  Quincke  Tube  - Continued 


First  Card 


COLUMNS  FORMAT 


NOTE:  If  there  are  no  holes,  only  the  first  card  is  required. 


i*.  :*  :•»  :i  ji'  a j*  jc  it  • 40  «i  «2  45  44  -it  «o  4 ■ j o jo  94  *.;•  n 9«  • vj  • 1 to  c c?  tl  r.«  | 


1 yi  ji  y/  :i  J4  J9  -.c  it  no  59  40  4i  4?  4j  . 


0 0 0 G 0:0  C 0 0 | 0 0 0 3 C 

1 I I I 1U  I I I '!  u II  1!  1 

1 I t 1 1 1 1 I I I II  I I 1 1 

2 2 7 7 22Z7T222? 

3 3 1 3 3 3 3 3 1 3 1 3 3 1 3 
444444444114441 

1 i 1 • i[»  1 1 t 1!  1 w 1 . 

5 5 b 5 b'b  b 5 5 5 5 5 5 5 


3 1 3 c c o c 0 0 o;o  0 o 3 cju  o c o go  o ||c  o 3 o o c;o  o o o o 3 o o o o;o  1 1 1 o 3 o o o o;o  0 0 0 0 3 0 D 0 0;0  0 0 3 3 


i mi  .i  • i:  h :i  . 

I I I I I I ! ! I I 
:727727222 
1333333333 
4 4 4 4 4 1 4 4 4 4 

» 1 « i ?j  *i  ::  }j  :•  :j 

S555bb5555 


» H ii  :i  y.  it  i j4  . i'  ,i  ii  4,  ii  4?  (j  t4  t*  n v » '.]  u v> 

1 1 1 1 1 1 1 1 1 pi  1 1 1 1 1 1 it  i|i  1 1 1 1 1 1 1 1 1 

i i 

222722222  7'2  22227222  2‘2  222722222 

I I 

1 3 1 3 3 3 3 3 3 3 | 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 

4 1 4 4 4 4 4 4 4 4 '4  4 4 4 4 I 4 4 4 4<4  4 4 4 4 4 4 4 4 4 

‘ .’I  l lijii  ! u t l<  ii  ,i  tl  ii  *!  il  n i . 1 1 1 1 ii  i»  u ii  ,i  -.1  :i  ii 

b b 5 b 5 b b b b 5 5 |b  5 5 5 b b 5 bib  5 5 b 5 b 5 5 5 b 


C G C G S S £ 6 6 ('c  6 6 E | i £ 5 G G'O  S 5 i f G C £ |c!i  S 0 E 6,5  6 G G GIG  G E £ C.£  6 & | 6 6 6 6 6 6 i G G 6 G 5 6 G 6 E’G  6 G 1 6 3 5 6 6 G 6 G 6 G 6 

UCCONMCLL  AUTOMATION  COMPAN  ■ [ [ 

7 2 1 1 i;i  I 7 2 I T 1 I T 7 I 7 I ; Jil  ? 2 1 T'T  I 7 I Til  I 7 I I'l  1 I I 1 I I J 1 T J I ? 7 7 I 7 T 2 7'l  J ) 7 T 1 T 7 7 I I 7 T 7 1 7 J 2 1 7'J  ) 7 1 J 

iii  i i i i 

'i'i  i i i i 

I I I ! | ! ! J 3 I 4 J H 1 1 1 l U I I t 3 : 1 3 4 S 6 3 3 4 I S 1 1 S I I II  . I S ||  3 ! 8 3 8 8 S 6 1 8 t ! 4 I J 8 8 5 1 8 I 8 6 5 S S S 8 1 8 ii  8 

; ; I I ; j I ; ' I i 

3 1 1 JU  3 5 3 9 3 4 3 4 1 i 9 9 13 ,3  3 3 9 s|j  3 1 3 9 13  4 1 0 VG  3 5 3 4 3 3 '■  jS  1 |j  9;9  9 ‘3  1 3d  3 ! 5 3 G G 9 5 U 3 I 3 j'l  3 9 J j'j  4 3 3 3 


,>i  H b)  i:  •)  (?  a (4  l.,U  t:  (i  (i  it  n 12  ii  i .»  i 'i  '1  ii 

; i ii  1 1 1 1 1 1 1 ; 1 1 1 1 1 1 1 1 1 1;  1 1 1 1 1 

2 2 2 2 3 2 2 7 2 2'2  2 2 2 2 2 7 1 2 J‘2  3 2 3 1 

I I 

( 3 3 3 3 3 3 3 3 3 1|  1 3 3 3 3 3 3 3 3;|  3 3 3 3 

4 4 4 4 4 4 4 4 4 i;4  4 4 4 4 4 4 4 4 4 4 | 4 4 4 

it  I!  il  II  it  II  I’  tl  II  tl ; it  I'  tl  fl  It  II  11  11  <1  Itj  't  t I II 

b 3 5 5 5 b 5 5 5 bjb  3 b 5 b b b b b 5,5  3 5 b 5 

I I 

6EGEC6EG6  C'E  6C|e  ioEECGCSGS 

I I 

f 1 1 1 1 1 1 1 1 ?;!  I 1 1 1 1 1 1 M'l  I 1 1 1 
|t  ? J J 8 8 8 8 8 Is  * 8 6 c : t S s|s  5 t 8 


» III  I •«  b I 


I * .4  I 1/  II  I I-.  V .r  3|  »»  * 
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2.3.8  Quincke  Tube  - Continued 


Second  Card 


COLUMNS  FORMAT 


DIMENSIONS 


11-20 


F10 . 0 Number  of  holes 

F10.0 Length  of  holes 


H.  n . 03d 

..  t 7 i I & r v >i  io  i.  a B « ii  i.  <t  iti  . . .i  t:  . J ; » : . : - *•  :*  :*  so  ;i  97  oi  »•  K.  36  it  it  ;».o  .<  .7  «j  .«  sc  .7  «j  49  w 51  17  w M s„  u m ’-9  to  t<  c:  C9  t. 

• t » • * *•  f • » i3  H U U u ij  il  if  it  r*  Ij  i < 14  i\  tt  ; ?■)  jo  v y:  ;j  J4  41  3b  yr  M yy  4j  «i  4 i 43  44  44  «6  4f  43  43  to  it  5.’  43  14  v.  M t?  tj  43  *o 

0 0 0 0 C;3  1 0 0 0 0 0 0 0 B O 1 3 0 S;3  0 0 0 0|0  C II  D oL  3 D 0 O'!  0 t C 0 0 0 0 0 C;0  0 0 G B 'o  0 0 0 O'O  0 0 0 0 jo  0 0 0 0,0  0 0 0 ojo  J 0 0 [TO  0 0 0 0 

1 l »*  t 1 it  i i;  ii  :«  i.jn  tin  •;  ;tjji  n n n ii|?k  i n :i  « r .1  :i  i»j:>  • ’»  !i  r n 1?  4)  1 4\,jt  u 11 it  m hJm  ti  4 •» hi  i?tl  11  o,ii  ft  ti  is  ;«  » u u ::  1:311 

M I 1 l!l  I I I 1 I I I I M 1 I I I I i ! I I'l  II  II  M I I I'l  I I I I I II  I l!l  I 11  I II  1 I I'l  I I I I I I II  I'l  1 I 1 I I I II  I'l  I I I I 


’i  :i  rjn  1?  4j  1 in, 


I I I I III!  I I I'l  I I I III  :!  I M I I I I'l  I I ! I'l  I I I III  I 1 I I 


2 2 2 2 2 2 2 2 2 2 '?  2 2 | 2 2 2 2 2 2 ? 2 2 2 2 2 ???  2 
1 • 

| j 3 3 3 3 n 33^  1 1 i n 3 ] no 3 3 j 3 3 3 3 3 3 

i I I 4 t I • 4 4 I'M  4 4 * 4 I 4 4 Lit  4 4 4 I 4 I 4 4 

* I » I • ? 1 14  » I II  '3  H 1 " .1  J>j’»  I'  .»  ’!  i.  I J’  .]  14  .) 

5 5 5 5 i i i'5  5 5 5 9 } 5 i 5 i'S  S i i i 5 1 i i i 


3 2 2 2 1 J 3 3 2 20  ? 2 2 2 2 2 3 J 2'2  2 2 2 JO  1 2 1 2'2  2 2 2 2 7 3 2 2 3'2  2 2 2 2 2 3 2 2 2 1 2 2 2 2 

; | | i i 

333333333  30  33333333  3*3  333  3 3 333  3'3  33333333  3'3  33333333  3'3  3333 

: I I i : | i 

niHonuiHtMuunMuniHiKnm'tniniumm 


j ' :i  *,  *i  u *i  t 
•i  i i 5 5 i 5 3 ‘j  i 


( i t C 6'5  G 6 6 C C 6 6 6 6 S 5 6 6 6lS  5 C 6 S S 6 S 6 CI3  C C 6 C 5 E C 6 G,6  £ C f 6 

I I MCDOW/CLL  AUT  OMATIC/N  COMI->AHY 

m 7 1;?  2 1 j 1 1 m i'l  n i i n i j I'l  1 1 / 1 1 / 1 ? jo  1 1 ; i n m 
• ; : : 

4 t 3 ! || ! ! ! ! I 6 I I ! | ! ! II 8 8 I I M I l It il I I LI  I I I 13  I I I I 

18  II  l|j  8 I 13  M 81  il)  5 3 9 9 9 9 8 3 9 0 9 3 9 1 3 9 13  l|l  9 8 I 8 ! J I 9 9 


i S 5 1 Mi  1 S i i'i  > i 5 5|S  5 5 5 J'i  5 i 5 i,5  5 5 i i'5  5 5 5 5 


C t C 5 G 5 6 6 6 6 6 C 6 C E S C G G SO  C C 3 E $ 6 C 5 C C 6 C 6 E 
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2.3.8  Quincke  Tube  - Continued 
Third  Card 


COLUMNS 

FORMAT 

DATA 

DIMENSIONS 

1-10 

F10.0 

Distance  to  first  hole 

INCHES 

11-20 

F10.0 

Distance  to  second  hole 

INCHES 

21-30 

F10.0 

Distance  to  third  hole 

INCHES 

31-40 

F10.0 

Distance  to  fourth  hole 

INCHES 

41-50 

F10.0 

Distance  to  fifth  hole 

INCHES 

51-60 

F10.0 

Distance  to  sixth  hole 

INCHES 

61-70 

F10.0 

Distance  to  seventh  hole 

INCHES 

71-80 

F10.0 

Distance  to  eighth  hole 

INCHES 
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2.3.8  Quincke  Tube  - Continued 


Fourth  Card 


COLUMNS 

FORMAT 

DATA 

DIMENSIONS 

1-10 

F10.0 

Diameter  of  first  hole 

INCHES 

11-20 

F10.0 

Diameter  of  second  hole 

INCHES 

21-30 

F10.0 

Diameter  of  third  hole 

INCHES 

31-40 

F10.0 

Diameter  of  fourth  hole 

INCHES 

41-50 

F10.0 

Diameter  of  fifth  hole 

INCHES 

51-60 

F10.0 

Diameter  of  sixth  hole 

INCHES 

61-70 

F10.0 

Diameter  of  seventh  hole 

INCHES 

71-80 

F10.0 

Diameter  of  eighth  hole 

INCHES 

NOTE:  If  required,  the  fourth  card  is  distances  for  holes  9-16,  with  the  same 

format  as  the  third  card.  If  number  of  holes  is  1-8,  the  fourth  card  is 
diameters  of  holes  1—8,  as  shown. 
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2.3.8  Quincke  Tube  - Cont'd 


Fifth  Card 

If  the  number  of  holes  is  1-8,  this  card  is  omitted. 

If  the  number  of  holes  is  9-16,  this  card  contains  the  diameter 
of  holes  1-8  as  shown  above. 


Sixth  Card 

If  the  number  of  holes  is  1-8,  this  card  is  omitted. 

If  the  number  of  holes  is  9-16,  this  card  contains  the  diameter 
of  holes  9-16. 
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The  hydraulic  pump,  i.e.  the  acoustic  noise  source,  is  always  treated 
as  the  first  element  in  the  circuit,  and  is  identified  as  a NTYPE  "9"  element. 


However,  the  program  is  capable  of  handling  different  pump  subtypes.  The  user 
must  select  the  pump  model  to  be  used,  and  input  data  accordingly,  as  described 
below.  Two  basic  pump  models  are  presently  available,  a complete  model  and  an 
empirical  model. 

The  complete  pump  model  (PUMP  subroutine)  is  based  on  actual  physical 
t dimensional  data  of  the  pump.  The  empirical  pump  model  is  based  on  experience 

gained  while  using  the  complete  model,  and  may  be  used  to  save  cost  or  until 
hardware  data  is  available  on  a new  pump.  It  may  also  be  used  to  model  other 
types  of  noise  sources  whose  characteristics  are  known.  Physical  data  for  a given 
pump  is  read  into  the  element  data  list  in  the  same  manner  as  for  the  other 
system  elements.  Required  input  data  for  the  complete  and  empirical  pump  models 
are  described  in  the  following  paragraphs  and  figures. 

The  complete  pump  model  has  three  subtypes  available  for  use.  KTYPE  21 
is  used  when  analyzing  the  pressure  side  of  the  system,  i.e.  no  return  side 
model.  KTYPE  22  is  used  for  pressure  side  analysis,  but  also  provides  a limited 
pump  inlet  analysis  sufficient  for  studying  pump  hanger  torque  characteristics 
without  the  need  for  a return  system  model.  KTYPE  23  is  used  when  analyzing 
both  the  pressure  and  return  sides  of  the  system. 

Data  for  the  pump  is  the  same  for  all  three  pump  subtypes.  However, 
the  system  data  must  contain  only  the  pressure  side  elements  for  KTYPE  21 
and  22. 

2. 3. 9.1  Complete  Pump  Model  - Input  data  for  the  complete  pump  model  requires 
three  data  cards  in  the  sequence  described  below.  Physical  data  are  input  for 
sixteen  variables  in  the  PUMP  subroutine.  Required  input  data  are  described 
in  the  following  tables  and/or  diagrams. 
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2. 3. 9.1  Complete  Pump  Model  - Cont  d 


COLUMNS  FORMAT 


DIMENSIONS 


11-20 


21-30 


31-40 


41-50 


51-60 


61-70 


71-80 


F10.0 


F10.0 


F10.0 


F10.0 


F10.0 


F10.0 


F10.0 


KTYPE 


NTYPE  = 9 


21  Pressure  Acoustics 


KTYPE  = 22  Pressure  Acoustics, 


Hanger  Torque 


KTYPE  = 23  Pressure  and  return 


acoustics,  hanger 


torque 


R1  = Cylinder  slot  radius 


SLOTW  = Cylinder  slot  width 


RV  = Cylinder  and  valve  plate 


slot  centerline  radius 


RB0RC  = Cylinder  centerline  radius 


DIAPIS  = Piston  diameter 


POVOL  = Oil  volume  between  piston 


at  mid-stroke  and  port  face 


R2  = Valve  plate  pressure  slot  radiu: 


INCHES 


INCHES 


INCHES 


INCHES 


INCHES 


INCHES 
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2. 3. 9.1  Complete  Pump  Model  - Cont'd 


Second  Card 


COLUMNS  FORMAT 


1-10  F10. 0 


11-20  F10.0 


21-30  F10. 0 


DIMENSIONS 


R4  = Valve  plate  suction  slot  radius  INCHES 


SWASH  = Maximum  swash  angle  DEGREES 


TLEAK  = Pump  internal  leakage  at 


input  steady  state 


pressure 


ANGCR  '=  Swash  plate  fixed  cross 


DEGREES 


swash  angle  plane) 


THPRS  = Valve  plate  pressure  slot 


start  angle 


THPRE  = Valve  plate  pressure  slot 


end  angle 


THSUCS  = Valve  plate  suction  slot 


start  angle 


THSUCE  = Valve  plate  suction  slot 


end  angle 


DEGREES 
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2. 3. 9.1  Complete  Pump  Model  - Cont'd 


Third  Card 


COLUMNS  FORMAT 


F10. 0 


Ll-20 


21-30 


F10.0 


F10.0 


31-40 


F10.0 


DATA 


LPRESS  = Pump  suction  port  (inlet) 


steady  state  pressure 


HOFF  = Swash  plate  centerline  offset 


DISAM  = Maximum  Swash  Plate  Actuator 


displacement 


ACTLEVO  = Swash  plate  actuator  lever 


41-50 


51-60 


61-70 


F10.0 


F10.0 


F10.0 


PIMASS  = Pumping 

piston  mass 

CPRESS  = Steady  state  case  pressure 

CSPRESS  = Case  to 

inlet  pressure  differ- 

ence  at 

zero  case  drain  flow 

DIACT  = Diameter 

of  swash  plate 

actuator 
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Mid  stroke  oil  volume  (POVOL)  is  the  volume  between  the  port 
face  and  the  bottom  of  the  hollow  piston. 


FIGURE  2 5 

SWASH  PLATE  PARAMETERS 


2. 3. 9. 2 Empirical  Pump  Model  (KTYPE  = 0) 


Observations  of  the  complete  pump  model  show  that  the  pump  flow  and  shunt 
impedance  can  be  approximated  using  values  of  the  variables  calculated  by  the 
complete  model.  Both  the  flow  and  the  shunt  impedance  are  kept  in  complex  form  to 
maintain  the  phase  relation  between  these  variables  and  pump  rotation,  as  predicted 
by  the  complete  model. 

Pump  Flow  (QO) 

An  oscillatory  flow  output,  QO,  of  the  "constant  current"  source,  is  assumed 
to  exist  within  the  pump,  and  is  represented  by  a complex  number.  The  fundamental 
component  of  the  pump  dynamic  flow  output  at  2500  rpm  is  quoted  in  cubic  inches  per 
second,  rms.  The  program  assumes  the  oscillatory  flow  to  be  proportional  to  pump 
speed,  and  calculates  the  input  flow  applicable  to  each  rpm  from  the  flow  at  2500 
rpm. 

If  the  flow  values  could  be  obtained  experimentally,  the  imaginary  part  of  the 
complex  flow  could  be  zero,  and  the  complex  flows  and  pressures  obtained  for 
points  elsewhere  in  the  system  would  be  phase  referenced  to  this  flow.  How- 
ever, if  the  flow  figures  are  obtained  from  a prior  program  run  in  which  the  complete 
pump  program  was  used,  both  real  and  imaginary  parts  will  result  as  this  flow, 
and  all  quantities  calculated  by  the  program  will  be  phase  referenced  to  valve 
port  opening. 

Pump  Shunt  Impedance  (20) 

An  impedance,  ZO  is  assumed  to  shunt  the  "constant  current"  flow  source 
at  2500  rpm,  and  is  represented  by  a complex  number  (psi  rms/cis  rms).  The  program 
currently  assumes  that  the  impedance  input  is  proportional  to  the  pump  rpm,  the 
fluid  viscosity,  and  the  fluid  mass  density.  The  shunt  impedance  for  each  pump 
rpm  is  calcaulated  from  the  shunt  impedance  at  2500  rpm.  This  relationship  is 
derived  empirically. 
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The  pump  output  pressure,  PO,  is  represented  by  a complex  number  (psi  rms) . 
This  input  provision  is  intended  to  be  used  when  test  results  from  the  circuit 
under  investigation  are  available.  Changes  must  be  made  to  the  program,  however, 
to  allow  this  input  to  be  used. 
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2. 3. 9. 2 Empirical  Pump  Model  - Cont'd 


-5 


6-10 


11-20 


21-30 


31-40 


41-50 


FORMAT 

DATA 

15 

NTYPE  = 9 

15 

KTYPE  = 0 

E10.3 

Pump  Flow  (QO)  - real  input 

E10.3 

Pumf  Flow  (QO)  - imaginary  part 

E10.3 

Pump  Shunt  Impedance  (ZO)  - real  part 

E10.3 

Pump  Shunt  Impedance  (ZO)  - imaginary 

DIMENSIONS 


51-60  E10.3 


61-70  F5 . 


Pump  Pressure  (PO)  - real  part 


Pump  Pressure  (PO)  - imaginary  part 
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2 . 4 Output  Requirements 


Input  data  cards  are  required  in  order  to  specify  the  desired  program 
output  results,  which  may  consist  of  five  categories:  flow  plots,  pressure 
plots,  impedance  plots,  acoustic  energy  density  plots,  and  acoustic  intensity 
plots,  vs.  pump  speed  (rpm) . At  least  one  card  is  required  in  each  category 
to  specify  the  total  number  of  plots  desired;  and  each  location  (element  number) 
in  the  circuit  for  which  a plot  is  to  be  generated.  If  no  plots  are  desired 
in  a category,  a blank  card  must  be  included,  in  the  proper  order,  for  that 
category.  All  cards  in  this  output  group  must  be  in  the  order  specified  below. 

Each  plot  card  contains  a series  of  integer  numbers.  The  first  gives  the  total 
number  of  plots  required  in  that  category.  The  remaining  numbers  given  the  identi- 
fications of  the  circuit  element  numbers  for  which  plots  are  to  be  produced.  For 
instance,  a flow  plot  specified  for  circuit  element  number  6 (line)  in  the  system 
described  in  Section  2.1,  produces  a flow  plot  at  the  junction  of  elements  6 and  5; 
i.e.  the  input  flow  (upstream  end)  to  the  specified  element  (6).  All  input  data 
numbers  occupy  five  columns  so  that  the  eighty  columns  of  a data  card  can  contain 
sixteen  numbers,  the  number  of  plots  plus  fifteen  element  numbers.  When  the  card  is 
read,  the  program  reads  the  total  number  of  plots  required  in  that  category  and  then 
looks  for  that  number  of  element  identifications.  If  the  number  is  greater  than  15 
the  program  starts  to  read  the  next  card.  For  example,  in  the  pressure  plot  cards 
shown  as  an  example,  since  the  total  number  of  pressure  plots  called  for  is  17,  the 
program  reads  the  first  card  and  then  completes  its  list  of  elements  by  reading  the 
numbers  it  finds  in  the  first  two  five-column  fields  on  the  second  card.  Example 
pressure  cards  are  shown  for  the  illustrated  F-15  PC-1  system  input  data. 


2.4  Output  Requirements  - Cont'd 


First  Card 


Pressure  Plots  - First  Card  (Typical) 
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2.4  Output  Requirements  - Cont'd 


Second  Card 


FORMAT 

DATA 

15 
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2.4.1  Pre-deflned  Output 

Any  of  the  normal  five  plot  categories  may  be  used  to  plot  other  variables 
computed  within  the  program.  This  is  useful  for  studying  the  behavior  of  certain 
variables  as  a function  of  pump  speed.  Normal  output  titles  are  retained  when 
making  special  plots. 

Eight  variables  associated  with  pump  performance  are  currently  pre-defined 
in  the  pump  subroutine.  Any  or  all  of  these  variables  will  be  plotted  if  they 
are  selected  in  the  normal  manner  as  described  in  2.4  above.  These  variables 
have  been  assigned  in  the  last  four  positions  of  the  flow  (Q)  and  pressure  (P) 
arrays.  The  position  and  description  of  the  pre-defined  variables  are  as 
follows . 


Assigned 

Name  and  Location 


Description 


Units 


Q(37) 

Q ( 38  ) 


Q ( 3 9) 


Q(40) 


P (37) 


P (38) 


P (39) 


P (40) 


Swash  plate  angle  DEGREES 

Pre-compressed  piston  pressure 

minus  steady  state  outlet  pressure  PSI 

(absolute  value) 

Piston  Pressure  at  start  of  decom-  PSI 

pression 

Piston  cavitation  at  end  of  decom-  INJ*1000 

pression 

Decompressed  piston  pressure  minus  PSI 

steady  state  'diet  pressure  (absolute 
valve) 

3 

Piston  cavitation  at  start  of  pre-  IN  *1000 

compression 

Piston  pressure  at  start  of  pre-  PSI 

compression 

Swash  plate  torque  due  to  piston  IN-LBS 

pressure  and  piston  acceleration  forces 


The  above  plot(s),  if  selected  for  output,  will  plot  regardless  of  the  number  of 
circuit  elements  in  the  model. 
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2 . 5 Program  Deck  and  Input  Data  Job  Set-Up 

The  program  deck  and  input  data  cards  are  assembled  as  shown  in  Figure  2-6. 
All  input  data  cards  must  be  in  the  order  shown.  Circuit  element  data  cards 
must  be  sequential  starting  with  element  number  1 and  ending  with  the  last, 
highest  number  terminated  element.  Job  control  card  content  and  arrangement 
varies  depending  on  the  particular  computer  facility  available  to  the  user. 


i 

l 


■ .. 
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FIGURE  2-6 

HSFR  PROGRAM  DECK  AND  INPUT  DATA  JOB  SETUP 
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3 . 0  PROGRAM  OUTPUT 


3 . 1 Computer  Output 

Figure  3-1  shows  example  outputs  printed  by  the  program  prior  to  the  output 
plots  for  the  illustrated  F-15  PC-1  system. 

The  program  automatically  prefaces  the  output  with  the  program  name,  date  of 
run,  run  title,  pump  speed  range  and  increment,  harmonic,  and  a list  of  the  fluid 
properties  calculated  by  the  program  for  the  specified  fluid,  temperature  and  pres- 
sure. Then,  a list  of  all  the  elements  used  in  the  run,  together  with  their  input 
data,  is  printed.  The  program  then  prints  the  total  number  of  flow,  pressure, 
impedance,  acoustic  energy  density,  and  acoustic  energy  intensity  plots  desired 
and,  in  each  case,  the  specified  circuit  element  numbers  for  which  plots  are 
desired . 

In  the  absence  of  special  instructions,  the  output  print  plots  are  printed 
next.  Extra  WRITE  statements  may  be  incorporated  for  printing  out  any  desired 
parameter  contained  in  the  program.  Extra  WRITE  statements  are  useful  to  investi- 
gate the  behavior  of  certain  variables. 

3 . 2 Output  Plots 

The  normal  computed  program  output  consists  of  plots  of  flow,  pressure, 
impedance,  and  acoustic  energy  density  and  intensity  versus  pump  rpm,  for  the 
selected  junctions  in  the  system  being  analyzed.  The  manner  in  which  the  user 
specifies  the  plots  required  is  described  in  Section  2.4.  The  functional  sub- 
routine GRAPH2  used  to  produce  the  plots  is  described  in  Volume  IV,  and  the 
restrictions  placed  by  GRAPH2  in  the  specified  pump  speeds  are  given  in  Section  2.2. 

3 . 3 Interpretation  of  Computer  Output  Plots 

A typical  computer  print-plot  output  shown  in  Figure  3-2  is  a plot  of  the 
first  harmonic  (fundamental)  frequency  amplitude  of  the  oscillatory  pressure,  versus 
pump  speed  (rpm) , for  element  number  7 in  the  illustrated  F-15  PC-1  system.  The 
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FIGURE  3-1  Typical  Pre-Plot  Computer  Output 
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pressure  amplitude  is  the  peak  value  of  the  pressure  sine  wave  for  the  fundamental 
frequency  only,  at  the  junction  of  elements  6 and  7.  This  plot  is  synonomous  with 
test  results  obtained  when  using  a spectrum  analyzer  storage  oscilloscope  to  record 
fundamental  frequency  pressure  amplitudes  at  that  junction  in  the  system  as  pump 
speed  is  varied  over  the  rpm  range  (photograph  in  Figure  3-2) . 

Pump  speeds  which  cause  maximum  amplitude  resonant  response  in  the  circuit 
are  identified  by  the  major  peak  amplitudes,  in  this  case  at  about  1450,  2850,  and 
3125  rpm.  Note  the  numerous  minor  resonant  responses  which  occur  in  a complex 
multi-branch  circuit  such  as  the  illustrated  F-15  PC-1  system.  For  a given  major 
resonant  speed  (rpm) , the  computed  peak  input  pressure  at  each  element  may  be 
plotted  (manually)  against  the  physical  location  of  each  element  in  the  circuit. 

The  resulting  standing  wave  pattern  is  shown  in  Figure  3-3,  which  compares  computed 


pressures  to  actual  test  measurements  made  on  the  F-15  iron  bird  PC-1  hydraulic 
system.  This  form  is  useful  in  analyzing  potential  system  modifications  which 
may  reduce  and/or  relocate  major  resonant  conditions. 

A similar  computer  output  plot  and  standing  wave  plot  are  shown  in  Figures 
3-4  and  3-5,  respectively,  for  a simple,  straight  test  circuit  modeled  and 
tested  with  the  F-4  aircraft  Abex  pump.  Note  the  lack  of  minor  resonant  responses 
in  the  simple  circuit  and  the  excellent  prediction  of  both  the  first  and  second  major 
resonant  speeds.  Fundamental  frequency  half  wave  amplitudes  predicted  by  the  earlier 
HSFR  program,  as  shown  in  Figures  3-2  through  3-5,  are  about  20%  lower  than  actual 


measured  vales.  Amplitudes  predicted  by  the  current  program  are  about  20%  higher 


than  measured  values.  This  is  primarily  due  to  changes  made  in  the  PUMP  subroutine. 

It  is  necessary  to  understand  the  behavior  of  acoustic  waves  in  a closed 
system  in  order  to  interpret  the  relationship  between  the  pressure  and  flow  at 
a particular  point  in  a system,  and  the  amplitudes  at  other  positions.  A 
typical  system  will  contain  a series  of  traveling  and  standing  waves  for  each 
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harmonic  frequency.  Appendix  A,  Volume  IV,  describes  how  these  waves  are 
distributed  along  a line,  and  the  relationship  between  pressure  and  flow,  and 
distance  along  the  line.  The  program  user  must  understand  this  part  of  the 
subject  to  preclude  drawing  erroneous  conclusions  from  the  results  of  a single 
point  pressure  measurement,  or  calculation. 

Acoustic  energy  analysis  may  also  be  performed  with  the  program.  Figure 
3-6  is  an  output  plot  of  the  acoustic,  energy  density  at  the  input  to  element 
number  6 in  the  F-15  PC-1  system.  The  energy  density  represents  the  total  of 
kinetic  (flow)  plus  potential  (pressure)  acoustic  energy.  In  this  case,  the 
total  energy  density  is  about  28  milliwatts  per  cubic  inch  at  the  low  resonance 
speed.  Figure  3-7  is  an  output  plot  of  the  acoustic  intensity  or  power  at  element 
number  6.  The  acoustic  intensity  is  analogous  to  active  AC  electrical  power 
flowing  past  the  junction  of  element  numbers  5 and  6.  The  active  power  is  product 
of  the  in-phase  component  of  the  oscillating  flow  and  pressure.  The  active 
acoustic  power  at  the  low  and  high  resonant  speeds  is  41  and  33  watts,  respectively. 
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FIGURE  3-2 

Computed  and  Measured  Fundamental  Frequency  Response,  F-15  PC-1  System 
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First  Harmonic  Peak 
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Spectrum  Analyzer 
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FIGURE  3-4 


Computed  and  Measured  Fundamental  Frequency  Response, 
Test  Stand  F-4  Abex  Pump 
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FIGURE  3-6 

EXAMPLE  OF  COMPETED  ACOUSTIC  ENERGY  DENSITY 
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FIGURE  3-7 

EXAMPLE  OF  COMPUTED  ACOUSTIC  ENERGY  INTENSITY  (POWER) 
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